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Abstract 

The issue of radiant spherical black holes being in stable thermal equilibrium 
with their radiation bath is reconsidered. Using a simple equilibrium statistical 
mechanical analysis incorporating Gaussian thermal fluctuations in a canonical en- 
t^j- ■ semble of isolated horizons, the heat capacity is shown to diverge at a critical value 

of the classical mass of the isolated horizon, given (in Planckian units) by the micro- 
canonical entropy calculated using Loop Quantum Gravity. The analysis reproduces 
<^ ■ the Hawking-Page phase transition discerned for anti-de Sitter black holes and gen- 

C***~ ■ eralizes it in the sense that nowhere is any classical metric made use of. 

o 

7 1 Introduction 

W). 

A fairly complete understanding of the entropy of isolated, macroscopic, generic four 
dimensional general relativistic black holes has been reached pQ within the framework 
of the isolated horizon [2] and its quantization in Loop Quantum Gravity [3]. Not only 
has the Bekenstein-Hawking area law been reproduced (albeit with a fit of the Barbero- 
Immirzi parameter) [3], and an infinite series of finite, unambiguous quantum corrections 
has been found [I], [5], with the leading correction, logarithmic in the horizon area, showing 
a degree of universality. 

Radiant black holes, in contrast, are not directly accessible by these methods. In- 
deed, isolated horizons, by virtue of their isolation, cannot radiate. However, a heuristic 
approach employing equilibrium statistical mechanics including Gaussian thermal fluctu- 
ations has been applied [6], [7]- [11], [12] to unravel the canonical (and grand canonical) 
entropy of radiating spherical black holes. Building on these results, this paper attempts 
to associate a phase structure with such black holes, and also to reveal the relevance of 
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the structure of quantum spacetime in arriving at this phase structure. Since no classical 
metrics are used anywhere in the analysis, in contrast to the incipient work of Hawking 
and Page [13] (on the thermal stability of spherical anti-de Sitter black holes), this work 
is to be thought of as a generalization of that work. The analysis makes crucial use of the 
existence of quantum corrections to the Bekenstein-Hawking area law for microcanonical 
entropy found within the Loop Quantum Gravity perspective [I]. 

The paper is organized as follows: in section II the quantum Hamiltonian constraint 
is shown to imply that canonical partition function can be expressed entirely in terms of 
a sum over states characterizing the event horizon (treated as a boundary of the black 
hole spacetime). In the continuum limit, the boundary partition function is expressed as 
an integral over horizon areas. The saddle point computation of the boundary partition 
function (with an isolated horizon of fixed the saddle point) is shown to lead 

to the same universal criterion for thermal stability obtained earlier by us [7] somewhat 
differently. In the next section, the phase structure of this thermodynamic system is 
investigated. It is shown that the transition from the stable to the unstable phase is a first 
order phase transition where the canonical heat capacity diverges. This phase transition 
can be thought of as a generalization of the first order phase transition between thermally 
stable and unstable black holes discerned by Hawking and Page [13] within a semiclassical 
analysis restricted to asymptotically anti-de Sitter spacetimes. In the concluding section, 
we point out certain exact properties of the boundary partition function, which has the 
potential to lead to a phase structure beyond the 'mean field theory' usually implied by 
the saddle point approximation. 



2 A Universal Stability Criterion 

To begin with, we adhere to the standard definition of a classical black hole spacetime as 
B = M. — J~(Z + ) where Ai is the entire spacetime and the J~ is the part of spacetime 
in the chronological past of future null infinity X + for asymptotically flat spacetimes. For 
asymptotically simple spacetimes like (anti)-de Sitter, the corresponding infinity is taken. 
Consider a canonical ensemble of such black hole spacetimes coupled to a radiation bath 
at a temperature The canonical partition function can be written as 

Z(0) = Trexp-(3H (1) 

where, the total Hamiltonian operator H can be decomposed as H — H v + Hf,. The trace 
is over composite (spacetime + matter) states which can be expanded over a generic basis 
set 

l*> = E^IVOItt), (2) 

u.b 

where \^ v ) represent composite states involving spin network-like states of bulk space 
within a Canonical Quantum Gravity scenario (like Loop Quantum Gravity for instance, 
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[E]) and bulk quantum matter states, while \\b) represent composite gravity- matter states 
associated with boundaries of spacetime. We assume that bulk states \ip v ) are annihilated 
by the bulk Hamiltonian 

H v \ip v ) = (3) 
as dictated by temporal diffeomorphism invariance. The trace then reduces to 

exp-pH bdy \xb) 
= Tr bdy exp -(3H bdy 

= Zbdy , (4) 

where, we have assumed that the composite matter-space bulk states are normalizable^ 
so that £„ |c^| 2 || \i[) v ) || 2 = £b G [0, oo). 

One concludes that bulk states have no direct contribution to the thermodynamics of 
gravitating systems with boundary; the latter is completely determined by states (degrees 
of freedom) on the boundary. One can think of this as a more general notion of hologra- 
phy than propounded earlier [T5] . The earlier idea essentially means that the boundary 
degrees of freedom contains the entire information of the bulk. But in diffeomorphism 
invariant systems, the bulk Hamiltonian must annihilate physical states, hence the entire 
information is what is encoded in the boundary degrees of freedom. Is there then no role 
of the bulk degrees of freedom ? The bulk states contribute as sources to the boundary 
dynamics. However, once this dynamics is determined by solving the boundary quantum 
equations of motion, this is sufficient to determine the thermodynamics of the system. 

Thus, in the particular case that the boundary in question is a black hole horizon, 
this is perhaps plausible evidence that its canonical entropy must be a function of some 
attribute of the horizon, like its area, although our approach does not actually prove 
that. To proceed further, we make one additional assumption : the eigenvalues of the 
Hamiltonian (which we call 'mass') associated with the boundary are a function of the 
eigenvalues of the area operator associated with the horizon, as determined within Loop 
Quantum Gravity. If g(M n ) is the degeneracy of the boundary state with eigenvalue M n 
, the boundary partition function can be written as 

Z{/3) = Y,9(M(A n ))e W -pM(A n ) , (5) 

n 

where, we focus on the large area eigenvalues of the boundary: A n = Airy/S'jlp ■ n for 
n> 1. As in [8], this enables us to write the partition function as an integral over the 
energy 

dM 



Z(J3) = J dMexp 



Smicro(M) — (5M — log 



dn 



(6) 



*It is not clear if this is a valid assumption in the case of Loop Quantum Gravity, see [14] . 
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we have used the standard relation between degeneracy of states and the microcanonical 
entropy, and the last term in the exponent is merely a Jacobian factor when one changes 
integration variables from n to M. Since dA(n)/dn = const., the partition function can 
be recast in terms of an integral over A 

POO 

Z(p) = / dA eiq>[S micro (A)-(3M(A)] (7) 
Jo 

The integral in (j7j) is performed using the saddle point approximation [8], with the 
saddle point chosen to be the area Am the equilibrium horizon which we assume to be an 
isolated horizon. Including the effect of Gaussian thermal fluctuations of the area around 
the equilibrium configuration, one obtains 



1/2 

Z IH ~ exp [S micro (A IH ) - (3M IH (A IH )\ 
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aa(Mih) 

where, the last square-root factor is the contribution of the Gaussian fluctuations around 
Ajh, with the subscripts implying derivatives with respect to area. Using the statistical 
mechanical definition of the canonical entropy, one gets [8J, [TT] 

S 'canon S m i cro — log A , (9) 



where, 



A = [M AA (A)S A (A) - M A (A)S AA (A)] , (10) 



where, k is a positive number independent of the area and we have dropped the subscript 
IH from M and A and micro from S which now represents the microcanonical entropy. 
The primes indicate derivatives with respect to the argument. It is clear that a sufficient 
condition for the isolated horizon being a point of stable thermal equilibrium is that 

A(A) > (11) 



which is precisely the condition found in |llj . 

The heat capacity can also be calculated and one obtains [TT] , using the microcanonical 
definition of temperature 

C = S\{A) A-\A). (12) 

This equation now leads to the same condition ( fTTT) as the necessary condition for ther- 
mal stability. Thus, spherical black holes are thermally stable iff ( TTTT) holds. It is not 
difficult to integrate this differential inequality with respect to area and obtain, with a 
few dimensional constants chosen appropriately for convenience, 

M(A) > S(A) . (13) 
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If we put back dimensional fundamental constants, this condition can be expressed as 
expressed as 



M(A) > 




(14) 



an inequality obtained earlier [TT] somewhat differently. 

3 A generalized Hawking-Page phase transition 

Clearly, C — > oo as A (A) — > thereby signifying a first order phase transition. It is easy 
to see that this phase transition happens for a critical area A c of the horizon, for which 



The canonical entropy, and hence the free energy have a discontinuity at this critical value 
of the area. This in turn can be used to define a critical mass M c = M(A C ). The variation 
of the heat capacity as a function of the equilibrium mass is depicted in Fig. 1. 



The phase transition thus occurs exactly when the mass of the equilibrium configu- 
ration, as a function of the classical area is identical (in Planckian units) with the mi- 
crocanonical entropy calculated from Loop Quantum Gravity and includes the infinite 
series of corrections over and above the area term [I]. In fact, the existence of at least 
the leading logarithmic correction with a negative coefficient turns out to be particularly 
important, for it ensures that Saa(A) > which is required for the existence of a phase 
boundary at the critical mass (see Fig. 2). 




(15) 





Fig. 1 Heat capacity as a function of equilibrium mass 
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Fig. 2 Phase diagram showing stable and unstable thermal equilibrium phases 

This phase boundary separates the phases corresponding to a black hole being in stable 
or unstable thermal equilibrium with respect to the radiation bath, given by the condition 
M(A) > S(A) or M(A) < S(A), as already discussed in [UJ. In terms of the effective 
microcanonical temperature, this corresponds respectively to T > I or T < I . 

In the foregoing analysis, no assumption has been made regarding the functional de- 
pendence of the equilibrium isolated horizon mass M(A); such a dependence of course 
originates from the behaviour of the classical metric corresponding to the black hole space- 
time. E.g., for the ordinary Schwarzschild black hole, it is known that M(A) ~ A 1 / 2 and 
therefore it does not satisfy our criterion above for thermal stability, as is well-known to 
be the case. The same is true for the Reisnner-Nordstrom black hole as can be easily 
seen by extracting the mass as a function of the horizon area for very large areas. We are 
of course assuming that the radiation from the black hole is composed of neutral, rather 
than charged, particles. In the latter event, the analysis will have to be done within a 
grand canonical ensemble [12] with perhaps similar end results. Thus asymptotically flat 
spherical black holes tend to have only an unstable thermal phase in which they radiate 
or accrete interminably. We are excluding here the possibility of the existence of yet an- 
other phase considered in [T3] where thermal stability of pure radiation against collapse 
into black holes has been considered. Our analysis is valid for spacetimes with an inner 
boundary and does not address the dynamical issue of how that boundary came about. 

Going over now to non-asymptotically flat geometries, the most important one is of 
course the adS variety, considered in detail in ref. [13J. From the adS-Schwarzschild 
metric, it is easy to read-off the mass-area relation 



where, I = (—A) 1 / 2 upto a numerical factor. The existence of two phases is obvious. If 
A is large, i.e., if the horizon area A > 4nl 2 substantially, to leading order in this area (in 
Planck units) M(A) ~ A 3/2 . Then, as is clear from Fie;. 2, one obtains stable thermal 
equilibrium. On the other hand, if A is smaller, such that A ~ Anl 2 , then M(A) ~ A 1 ! 2 




(16) 
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once again. One is then clearly in an unstable thermal equilibrium phase, since this 
is precisely the mass-area relation for an ordinary (asymptotically flat) Schwarzschild 
geometry, which corresponds to negative specific heat. Thus, the two phases precisely 
correspond to the two black hole phases discerned in ref. [13] for such a spacetime. Yet, 
nowhere in our analysis have we made any reference to classical spacetimes. 

What about the point of transition ? Our analysis, depicted in Fig. 1 would lead to the 
conclusion that the transition is first order at a critical mass equal to the microcanonical 
entropy as a function of the area. It is useful to compare this with the Hawking-Page 
transition from a stable adS-Schwarzschild black hole (for A » Anl 2 ) to an unstable 
one (for A ~ Anl 2 ). One can determine the precise point of transition in terms of A by 
referring to [6]; using eq. (28) of ref. [6], the heat capacity for an adS-Schwarzschild in 
four dimensions is given by 



C 



A 



M. + 1 

4tt/ 2 Z 

6A _ i 

iirl 2 L 



(17) 



The critical area is clearly A c = ^-l 2 - This is certainly consistent with that obtained using 
our formula ([TBI) and the adS-Schwarzschild mass-area relation (fTEj) . Similar conclusions 
emerge for the generic ads-Reissner-Nordstrom solution as well, of course restricted to 
the canonical ensemble for a neutral radiation bath. Once again, no recourse is anywhere 
taken to classical spacetime metrics in our analysis. 



4 Beyond the saddle-point approximation 

The saddle point approximation used in the analysis so far is similar to the approximation 
used in the incipient work of Hawking and Page to evaluate their partition function 
expressed as a functional integral over Euclidean metrics. The conclusions we have drawn 
therefore may be compared with theirs based upon their choice of the saddle point in the 
Euclidean functional integral to be the classical metric itself, together with appropriate 
boundary contributions. The improvement in our approach is that we have made no 
explicit reference to the classical metric, and hence have what may be thought of as a 
generalization of that earlier work. However, the use of the saddle point approximation, 
which may be thought of as 'a mean field theory' in the sense of many body theory, is 
often not considered reliable, since during a real phase transition, 'collective' effects may 
actually be quite large and thus lead to a rather different phase behaviour. If so, both 
the original Hawking-Page phase transition and our generalization of it will be subject to 
non-trivial corrections beyond the saddle point approximation. We have not made any 
substantive progress towards delineating the precise nature of these corrections within our 
approach. However, from a consideration of the starting formula (j7j) for the boundary 
partition function expressed as an integral over the horizon area, some facets of a phase 
structure appear to emerge without a detailed analysis. 
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Since both the mass and the microcanonical entropy are positive and monotonically 
increasing functions of the horizon area, it is obvious that the partition function in (JTj) is 
finite provided 

(5M(A) > S{A) (18) 

or, if 

" > s S <19) 

where the microcanonical entropy and mass are to be evaluated at equilibrium where one 
can assume that they refer to an isolated horizon for which both quantities are well- 
defined. 

Similarly, the thermally averaged mass, defined to be the equilibrium mass, is finite 
provided 

i.e., at a higher critical inverse temperature of the heat bath. Finally, the heat capacity 
is finite provided 

^ > fe ^ w+w* (21) 

an even higher inverse temperature. In this range, the heat capacity can be expressed as 

C = (M 2 )p - (M) 2 = AM 2 , (22) 

where, AM 2 is the mean squared fluctuation of mass, an obviously positive quantity, 
corresponding to a stable thermal equilibrium phase. However, for (3 < (3 C 2, the heat 
capacity appears to diverge, in general, thus strongly suggestive of a first order phase 
transition. 

The consistency of these general behaviour patterns with those obtained in the earlier 
sections through the saddle point approximation is not obvious. We hope to report on 
this and other aspects of these exact formulae elsewhere. 
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